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ABSTRACT 

We study numerically the interaction between an intense circularly polarized 
laser field and an electron moving in a potential which has a discrete cylin- 
drical symmetry with respect to the laser pulse propagation direction. This 
setup serves as a simple model, e.g., for benzene and other aromatic com- 
pounds. From general symmetry considerations, within a Floquet approach, 
selection rules for the harmonic generation [O. Alon et al, Phys. Rev. Lett. 
80, 3743 (1998)] have been derived recently. Instead, the results we present in 
this paper have been obtained solving the time-dependent Schrodinger equa- 
tion ab initio for realistic pulse shapes. We find a rich structure which is not 
always dominated by the laser harmonics. 

PACS numbers: 42.50.Hz, 33.80.Wz, 31.15.Ar, 42.65.Ky 
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I. INTRODUCTION 



The generation of harmonics (HG) through the interaction of atoms with intense laser fields 
is a topic which has been broadly studied from both a theoretical and an experimental point 
of view. The big interest in the HG is due to the possible use as a source of short-wavelength 
radiation. In fact, through the harmonic emission it is possible to generate coherent XUV 
radiation using table-top lasers. Recently, harmonics of wavelength as short as 67 A have 
been reported M. The harmonics spectra obtained from a single atom in a monochromatic 
laser field present some common and well known features: (i) only linearly polarized odd 
harmonics are generated f2j, (ii) the spectrum has a plateau structure , (iii) the plateau is 
extended up to a cut-off that is located around I P + 3.17U P , where I p is the ionization energy 
of the atom and U p is the ponderomotive energy ||. The presence of only odd harmonics is 
due to symmetry reasons (a more detailed argument will be discussed in the next Section) 
and the location of the cut-off can be explained, at least qualitatively, with the so-called 
"simple man's theory" ||. The interaction of a single atom with two circularly polarized 
lasers of frequencies u and 2 uj has been investigated recently || ; it has been found that while 
the harmonic orders 3 n — 1 and 3 n + 1 are allowed, the harmonic order 3 n is forbidden 
by the selection rules of the dipole emission, where n is any positive integer. These results 
agreed with a previous experiment ||. 

More recently, the generation of harmonics in more complex systems than the single atom has 
become a strongly addressed topic. A model for harmonic emission from atomic clusters has 
been proposed 0. The harmonics generated by electrons moving in the periodic potential 
of a crystal have been investigated also |||| . 

In this work we want to study the generation of harmonics in ring-shaped molecules, like 
benzene and other aromatic compounds. This kind of molecules exhibits an invariance under 
a rotation of a certain angle around the axis which is orthogonal to the molecule plane and 
goes through its center. In this case the potential is periodic in the azimuthal direction. 
The HG from ring-like molecules interacting with a circularly polarized field presents many 
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different features with respect to the single-atom case in a linearly polarized field: (i) within 
the same harmonic range fewer lines are emitted and the higher is the number of atoms 
in the molecule, the lower is the number of emitted lines, (ii) odd and even harmonics are 
equally possible, (iii) the harmonics are alternately left or right circularly polarized. In our 
opinion, all these peculiar properties make this topic challenging and worth to be studied in 
detail. 

The paper is organized as following: in Section II we summarize the derivation of the selection 
rules for the ring-shaped molecules obtained by Alon et al. [0 . In Section III the numerical 



model used in our simulations is presented and discussed. In Section IV we describe the 
interaction between the ring molecule and the laser field. In Section V we show the results 
obtained for different intensities and frequencies together with a broad discussion. Finally, 
in Section VI we give a summary and an outlook. Atomic units (a.u.) are used throughout 
the paper. 

II. SELECTION RULES 

In the case of atoms or molecules which are shone by a laser field of frequency u, the 
Hamiltonian is periodic in time with a period r = 2tt/uj: H(t + r) = H(t). The time- 
dependent Schrodinger equation (TDSE) for such a system can be written as 



* w (?,t) = 0. (1) 

The solutions are of the form (Floquet theorem): 

# w (f, t) = $(f, t)e~ iWt with $(r, t + t) = $(f, t), (2) 

where W is the quasi-energy and $(r, t) is a square integrable function. 
Because the set of all functions which are square-integrable in a certain interval and have a 
finite norm over a cycle, forms a composite Hilbert space, we can apply the extended Hilbert 
space formalism. The probability to get the nth harmonic from a system in a state \P}v(r, t) 
is gO | 



(n) 4 



(($(f,t)|/ie- m ^|$(f,t))) 



(3) 



where fi is the dipole operator and the double bracket stands for the integration over space 
and time. 

In the case of an atom (in the dipole approximation) the Hamiltonian Tt(t) = \H(t) — 
is invariant under the so-called second order dynamical symmetry operator (DSO) O], 

V 2 = {f -> -f,t ^t + n/u). (4) 

Therefore, the states |$)) are simultaneous eigenfunctions of 'Hit) and the second order DSO 
with eigenvalues ±1. The nth harmonic is therefore emitted only if 

(($|£ e -***|$)) = ((V2^\V 2 fte- iruJt V^\V 2 ^)) ^ 0, (5) 

leading to 

fi(r)e- inut = fi(-f)e- inuj{t+7T/uj) (6) 

that is fulfilled only with odd n's. Instead, where the Hamiltonian T~t(t) is invariant under a 
rotation around an iV-fold symmetry axis V 2 can be replaced by JIT] 



7» w =(v V ,+ — ,*-»* + — ) , (7) 

where ip is the angular coordinate around the symmetry axis. With an algebra similar to 
the single atom case we derive 

from which n = kN ± 1, k G J\f follows. That means that the higher is the symmetry order 
N the less are the generated harmonics within a fixed frequency interval. In the limit of a 
continuous symmetry a circularly polarized laser does not generate any harmonics. The 
two harmonics of each couple have opposite polarization, clockwise and anticlockwise QT2 . 
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III. NUMERICAL MODEL 



In order to keep the numerical effort manageable we restrict ourselves to a two-dimensional 
(2D) model where the molecule plane and the rotating electric field are properly represented. 
We study a single active electron in a ring-shaped potential of N ions. Different kinds 
of "smoothed" potentials can be used for this purpose [I3~|. The potential used in our 
simulations reads 

V(p,<p) = A (a cos (N ip) + 2 -a) (9) 

\/(p-r ) 2 + f3 

where r is the radius of the molecule, and p and (p are the polar coordinates, ft is the 
parameter which gives the "degree of smoothness" of the potential and determines the width 
of the ground state along the ring, a moves the local maxima of the potential keeping the 
minima constant (the parameter a was introduced in order to avoid the presence of points 
where V — for finite p, because this could generate non-physical ionization). Finally, A 
is the "strength" of the potential. For our simulations we chose a = 0.99 and (3 = 0.38. 
Once the values of a and (3, which lead to reasonable model properties, have been found 
and fixed, we varied A for choosing the ionization energy of the molecule. The potential has 
iV oscillations in the azimuthal direction, each minimum representing the location of one of 
the N ions. The potential goes to zero for p — > oo. 

It has been tested that the generation of harmonics is very weakly dependent on the fine 
details of the atomic potential. Instead, it is strongly dependent on its geometry and symme- 
tries. It is therefore worth to look for a model potential that, keeping the proper symmetry, 
can be quite easily numerically optimized, i.e., requiring as few as possible grid points. In 
order to achieve this minimization of the number of grid points and to not break any physi- 
cal symmetry we used a polar grid. Moreover, we paid attention to use always a number of 
points in the azimuthal direction that is an integer multiple of N. We will take N = Q, like 
in benzene, and therefore the potential will exhibit a Cq symmetry around the orthogonal 
axis. 
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For a good understanding of the harmonic spectra it is essential to study the level scheme 
in detail. Therefore, in order to characterize our model we have calculated the energy of 
the first six states for different potential strengths A in the interval between 0.2 and 1.6. 
In this interval the energy of the ground state of the molecule decreases from —0.14 down 
to —2.44. The first, the third and the fourth states are non-degenerate, the others have 
a double degeneracy. In Fig. 1 the energetic behavior of those states versus A is shown. 
Contour plots of the six states for an intermediate value of A = 0.80 are shown in Fig. 2. 
The pattern shape and the symmetries do not change with A, but the "average radius", 
i.e., the spatial extension, does. Clearly, for low values of A the probability density is more 
loosely bound than for high A. In particular this is true for the upper states. For the sake 
of easy reference later on, we name the six states 0, l a , 1&, 2 a , 2 b , 3, 4, 5 a and 5^. The 
subscripts are used just to distinguish among the degenerate states. The non-degenerate 
states 0, 3 and 4 have the full G% symmetry. 

As one can infer from Fig. 1 the six states can be divided in two branches: the first one, 
containing the first four states, — 3, decreases very fast as A increases, and a second one, 
containing the two upper states, 4 — 5, which decreases much more slowly and lies even 
in the continuum for A < 0.4. As a result, for an increasing A an increasing gap between 
the two branches appears. To obtain the ionization potential of real benzene (—0.34 a.u.) 
we have to choose A = 0.375. Surprisingly, at that position the level scheme of our simple 
model resembles the molecular orbital (MO) scheme of the real benzene very well [f[4| |. In 
particular, the states 4 and 5 are still in the continuum so that only the four states — 3, 
possessing the same degeneracies as the MOs, are bound. A magnification of the region 
around A = 0.375 is also shown in Fig. 1. Another parameter that will play a role in the HG 
is the level spacing Q between the ground state and the first excited state. For an increasing 
A, i.e., a decreasing ground state energy, fl decreases. 
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IV. MOLECULE-FIELD INTERACTION 



In dipole approximation the time-dependent Schrodinger equation for a single electron in a 
laser field E(t) and under the influence of an effective ionic potential V(f) is given in length 
gauge by 

i-*(r,t) = (--V 2 +V(r) + E{t) ■ r)*(r,t). (10) 

In our case the dipole approximation is excellent because the molecule has a size much 
smaller than the wavelength of the laser field. We used a circularly polarized laser field that 
in cartesian coordinates is described by 

E(t) = — j=-(cos(ujt) e x + siD.(ut) e y ), (11) 

where £(t) is a slowly varying envelope with amplitude £, and u is the laser frequency. In 
polar coordinates we obtain the TDSE 

9 T , \ ( 1 9 f d\ Id 2 rr , . S(t) . \ , . , 

l di^ ^ t) ={-2 pd - p {%)-2i 2 W 2+ K(P ' V) + 7f P C ° S( ^ " Ut) ) ^ *> t] • (12) 

This TDSE can be solved ab initio on a PC. We did this by propagating the wavefunction in 
time with a Crank-Nicholson approximant to the propagator U(t + At, t) = exp[—iAtH(t + 
At/2)} where H(t) is the explicitly time- dependent Hamiltonian corresponding to the TDSE 
(|l2l). Our algorithm is fourth order in the grid spacings Ap, A(f and second order in the 
time step At. The boundary condition is \I/(p, 0,t) = ^(p, 2tt, t) for all p and t. Probability 
density which approaches the grid boundary in p-direction is removed by an imaginary 
potential. 



V. RESULTS 

Here we discuss the results obtained from our 2D simulations and we compare them with 
previous results from a one- dimensional model (ID) presented elsewhere Our studies 
were mainly focused on the structure of the harmonic spectrum for different values of the 



ionization energy. In general, our findings show that, together with the harmonics we ex- 
pected from the selection rules, other lines are present and their location can be, in most of 
the cases, explained with the help of the level scheme. In particular, we observed that for 
higher A the gap between the two branches of Fig. 1 plays an important role. For each A var- 
ious simulations with pulses of the same frequency and length but different intensities were 
performed. We used sine-square pulses of 30 cycles duration and a frequency uj = 0.0942, 
unless noted otherwise. In order to better understand the additional lines which appear 
besides the expected harmonics it is useful to study the low intensity regime first. 



A. Low Fields 

In the case of a single atom, when the intensity of the field is not high enough for generat- 
ing harmonics efficiently, the Fourier transform of the dipole shows only the fundamental. 
Clearly, the threshold of the field strength to observe any harmonics depends on the ioniza- 
tion potential. In Fig. 3 two dipole spectra, for different A, obtained from the interaction 
of the ring molecule with a low field pulse are shown. For the dipole emission spectrum of 
Fig. 3a an electric field amplitude £ = 0.005 was used and A = 1.6. In this case apart from 
the fundamental and the weak fifth harmonic, four other lines are present. Two lines are 
located at 14.1a; and 16.0 uj, respectively. We refer to them as A a and A 6 hereafter. These 
two lines correspond to two resonances between the states 0-^5 and 3^4 and therefore 
they move towards the red if A decreases. This is confirmed in Fig. 3b, where A = 1.4 and 
£ = 0.005. In Fig. 3b the 5th harmonic is more pronounced. This is due to the fact that 
when the ionization potential is lower the generation of harmonics requires weaker fields. A 
subfundamental line at Q, corresponding to a transition from the first excited state 1 to the 
ground state is also present in both cases. We name this line T. Another line corresponding 
to the transition 2uj — Q is visible as a right satellite to the fundamental. Looking closer, 
one observes that, with respect to the positions one would expect from the unperturbed 
level scheme, A a and A b are blue-shifted, whereas T is not. This opposite shift can be 
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explained by the dynamical Stark effect. The ground state remains almost unaffected and 
the lower-lying states move slightly. The higher-lying states instead, experience a relatively 
strong shift. Therefore, with increasing laser intensity the gap between the two branches 
in Fig. 1 increases, leading to a blue shift of A a and A b . The three lines A a , A b and T are 
the first lines to appear in the low field regime. The subfundamental line can be considered 
as characteristic for the low field regime. Indeed, with increasing field strength, it moves 
towards the red and finally becomes very difficult to be resolved. 

B. High Fields 

With increasing laser intensity the actual harmonic spectrum develops. In Fig. 4 we show 
HG spectra for four different electric field amplitudes and A = 1.6. In Fig. 4a £ = 0.02 
a.u. and the situation is, at first sight, quite similar to that one of Fig. 3a: A a and A b are 
present, and the first allowed harmonic is there, i.e. the 5th. Increasing the field to £ — 0.12 
a.u., Fig. 4b, also the 7th and the 11th appear, and the two lines A a and A& are still present 
but located within a broader structure. This last phenomenon can be explained taking into 
account at least three different effects: (i) in general, the width of A a and A b increases for 
higher fields, (ii) other channels, i.e. other resonance lines between the two branches, can 
be opened if the field becomes intense, and (iii) through a removal of the degeneracy of the 
states by the electric field more possible resonance lines are obtained. These factors generate 
a kind of broad "hill" in the harmonic spectrum, the position of which is a function of A. 
Moreover, in Fig. 4b other satellite lines around the expected harmonics are present. Those 
lines are decays from virtual states to real states. In Fig. 4c, t = 0.21 a.u., more couples of 
harmonics are present and the hill is reduced to a light background modulation of the main 
HG spectrum. This becomes even clearer in Fig. 4d where harmonics up to the 49th are 
observed and the 17th is located just on the hill. We have also studied in more detail how 
the strength of the harmonics increases (or decreases) in function of the electric field. The 
results for harmonics up to the 31st are shown in Fig. 5. For each couple of harmonics there 
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is a minimum field threshold, below which the harmonic lines cannot be picked out from the 
background. All the harmonic strengths of Fig. 5 are normalized to the fundamental. What 
is worth to stress is that for the 5th harmonic we can have an efficiency up to 13%. 
Repeating the sequence of simulations shown in Fig. 4 for a lower A gives results that are 
quite similar to those shown in Fig. 4 as long as the gap between the two branches is quite 
large. When the gap becomes of the order of about ten photons the two lines A a and At play 
a role that is less important. This is due to the fact that, as the gap is smaller, the two lines 
are expected to be located in a low frequency region and therefore they are easily hidden 
by the background of the main HG spectrum that in the low frequency region is higher. 
Furthermore, as the intensity of the harmonic lines is strongly enhanced with increasing 
field, the strength of A a and A^ is not. Also the extension of the harmonic spectrum is 
dependent on A, for lower A less harmonics can be generated (for a fixed frequency). 
As already mentioned, Q is the distance between and 1. For high A (or high u) we have 
uj > Q but for decreasing A (or decreasing u) , Q approaches the laser frequency and overtakes 
it. The ratio of u and Q strongly affects the harmonic emission by the ring molecule. In 
particular, we observed that when u ~ Q a very complex spectrum is generated and together 
with the expected harmonics many other lines of similar intensities are present. This effect 



can be seen in Fig. 6. A very similar behavior was also observed in the ID model In 
particular, the shape of the additional satellite structures around each allowed harmonic 
are in the two cases alike. It seems that in this resonant case the system is not in a single 



non-degenerate Floquet state as assumed for the proof of the selection rules ||11|| . In that 
derivation the pulse was assumed as infinite. Therefore, the particular behavior of the dipole 
emission could be also due to a pulse shape effect. However, pulse shape effects should be 
not dependent on the frequency, and in fact, we have additional lines for all the frequencies, 
but those lines play always a minor role with respect to the expected harmonics. Instead, 
when the laser frequency becomes nearly resonant many new strong lines appear. 
Keeping the same pulse parameters and decreasing the parameter A the ionization increases. 
When we want to study a model that is closer to the real benzene, we have to take A = 0.375. 
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With this condition the physical scheme is very different from those of the cases previously 
discussed; the ionization energy is reduced from 2.44 to 0.34 and, as we already mentioned, 
only four states are present and the gap between level 3 and 4 does not exist at all. Under 
these conditions a pulse with the frequency uj = 0.0942, which we used so far, leads to 
emission spectra without any harmonic structure. This is mainly for the reason that the 
frequency is relatively high with respect to the ionization potential and therefore just a very 
few photons are sufficient for reaching the continuum. Therefore, unless the field is very 
low the ionization would prevail at soon. Also making a comparison with the rule for the 
cut-off position used in the atomic case (I p + 3.17U P ) we should not expect any harmonics 
due to the low value of U p . Therefore, we made a series of simulations with the benzene 
model but using a lower frequency, u = 0.0314. With this frequency eleven photons are 
required for reaching the continuum, i.e., the molecule can be ionized only with a high-order 
multiphoton process. A spectrum obtained with this low frequency and £ = 0.044 is shown 
in Fig. 7. Like in the highly-bound high frequency case the emission spectrum exhibits the 
harmonics allowed by the selection rules. The efficiency of the harmonics in Fig. 7 is not as 
high as the one of the harmonics of Fig. 4. Another line, located around 3 u is also present 
in the spectrum of Fig. 7. This line is a resonance between the states and 1 (for A = 0.37, 
Q = 0.0914 ph 3u ). It is interesting to note that in this case of a weakly bound electron the 
results from the 2D simulations are different with respect to those from the ID simulations 



15]. In the latter case no harmonic structure was observed for to < f2. This, in our opinion, 



could be due to the reason that in the ID model the level scheme is qualitatively different. 
In particular, there is no continuum in the ID case. 



VI. CONCLUSIONS 

In this work we have studied the harmonic emission in a ring molecule. We have shown 
that when a ring molecule interacts with a laser pulse, together with the series of harmonics 
predicted by the selection rules, other lines are present. Under certain conditions the strength 
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of these lines can be comparable with that one of the harmonics. Our HG spectra present 



a structure and a complexity that is absent in the numerical results shown in ||11|| . This is 
due to the reason that while there a ID Floquet simulation was performed in our studies a 
realistic pulse (i.e., a finite pulse with a certain envelope) and a 2D model with ionization 
included were used, and the TDSE was solved ah initio. 



What is worth to note is the scaling of the TDSE (fl2|) with respect to the size of the 



molecule. If one scales the molecule radius like p' = ap, the TDSE flT2| ) remains invariant 
if t' = a 2 t, V = V/a 2 , £' = S/a 3 , u' = uo/a 2 are chosen. Therefore our results for high A 
can reproduce the results that would have been obtained for a bigger molecule with a lower 
ionization potential interacting with a field of lower frequency. Moreover, one can think 
about the different cases of Fig. 1 as the level scheme of positively charged molecules as well. 
To our knowledge, so far harmonics in ring-shaped molecules have been investigated exper- 



imentally only in gaseous samples fl6| . Studying a gaseous sample is very different with 
respect to what we did in our simulations. Because the molecules in a gas have a random 
orientation it is not possible to apply the symmetry properties discussed. If the propagation 
direction of the circularly polarized field is not orthogonal to the plane of the molecule, 
the molecule "sees" a field that is elliptically polarized. This breaks the discrete rotational 
symmetry and other harmonics become possible. Nonetheless, even if it is not possible to 
make a direct comparison it is useful to note that the results presented in [|l6j confirm that 
ring molecules, like benzene and cyclohexane, can tolerate short pulses of high intensity and 
phenomena like fragmentation and Coulomb explosion do not play a big role. In order to re- 
produce in a real experiment the results we presented, it is fundamental to prepare a sample 
where most of the molecules lie in the same plane. This could be done with some orientation 
techniques or, considering the particular shape of the organic molecules, preparing a very 
thin layer. 

We described and discussed results for a molecules with N = 6, but cases with higher N 
are as well possible. Moreover, the higher is iV the higher is the frequency of the first 
generated harmonic; in the limit of very high N, harmonics of very short wavelength could 
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be generated. 

In this work we a took into account a single active electron, but also including correlation 



through both, a full description [T1J or an approximate treatment like time- dependent density 
functional theory, would not change the symmetry properties of the Hamiltonian describing 
the system. Therefore, the selection rules would apply as well. 

More complex molecules which produces similar selection rules [17]] are the nanotubes Jl9| . 



They can be very long in the longitudinal direction and exhibit a discretized cylindrical 
symmetry. A semiclassical approach to harmonic generation from nanotubes was also in- 



vestigated [18| . Unfortunately, because of their size, the dipole approximation would not be 
accurate enough. Therefore, an ab initio numerical simulation in 3D would be at the limit, 
or probably beyond, the calculation capabilities of even the fastest computers now available. 
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FIGURES 



FIG. 1. Energetic behavior of the first six states versus A. For lower A some of the states 
belong to the continuum; the magnification shows the region around A = 0.375 where four 
states are bound. For higher A the set of states is split in two branches and a gap between 
those branches appears. 

Fig. 2. Contour plots of the first six states. For each double degenerate state, two linearly 
independent states are shown. The non-degenerate states present a full C% symmetry. 

Fig. 3. Emission spectrum in the low field regime. In (a) A = 1.6 and £ = 0.005, in (b) 
A = 1.4 and £ = 0.005. The lines A a , A b and T are present in both pictures. When A 
decreases the two lines A a and A b are red-shifted. This can be observed comparing (b) with 
(a). 

Fig. 4. Evolution of the harmonic spectrum with increasing field for A = 1.6. In (a) 
£ = 0.02, in (b) £ = 0.12, in (c) £ = 0.21 and in (d) £ = 0.27. When the field increases 
additional satellite lines appear with the low order harmonics. The highest allowed resolved 
harmonic is the 49th. 

Fig. 5. Strength of each harmonic line versus the electric field £. It is worth to note that 
the first allowed harmonic, i.e., the 5th, can reach an efficiency up to 13%. In the low field 
region the strength of the the 13th and 17th harmonic is particularly high, this is due to the 
presence of the "hill". The value of A is constant, A = 1.6. 

Fig. 6. Emission spectrum for A = 1.00 and u = 0.0942. The value of Q approaches 
the laser frequency u and the spectrum exhibits a complex structure with many additional 
strong lines. 
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Fig. 7. Emission spectrum for A = 0.375 and 00 = 0.0314. Also in this case, where the 
ionization potential corresponds to that one of the benzene molecule, the spectrum exhibits 
the same structure as in the highly-bound high frequency case. The line around 3 u> is given 
by a decay from the state 1 to the state 0. 
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Fig. 2: F. Ceccherini and D.Bauer, "Harmonic generation in ..." 




Fig. 3: F. Ceccherini and D.Bauer, "Harmonic generation in ..." 




Fig. 4: F. Ceccherini and D.Bauer, "Harmonic generation in ..." 
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Fig. 5: F. Ceccherini and D.Bauer, "Harmonic generation in ..." 




Fig. 6: F. Ceccherini and D.Bauer, "Harmonic generation in ..." 




Fig. 7: F. Ceccherini and D.Bauer, "Harmonic generation in ..." 



